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^ ^ A method is developed to design a fixed-parameter 

D ^ compensator for a linear, time-invariant, SISO plant 
model characterized by significant structured, as well 
S as unstructured, uncertainty. The controller minimizes 
r? the H°° norm of the worst-case sensitivity function 
over the operating band and the resulting feedback 
( system exhibits robust stability ml robust 
performance. It is conjectured that such a robust 
nonadaptive control design technique can be used 
on-line in an adaptive control system. 

1. Introduction 

The mathematical description of a physical plant 
(i.e. the nominal model) is always characterized by 
uncertainty, or modeling error. The plant/model 
mismatch caused by neglecting high-frequency 
phenomena (i.e. unmodeled dynamics) is known as 
unstructured uncertainty because only a frequency 
dependent magnitude bound on the error is available. 
Differences between the actual values and the nominal 
values of the parameters in the finite-dimensional, 
low-frequency model are the source of structured 
uncertainty. The goal of the robust design method 
presented here is to maintain closed-loop stability and. 
performance in the presence of both types of 
uncertainty. 

The most popular modern design methodologies 
(LQG/LTR, H 00 ) only deal with plant models which 
contain unstructured uncertainty [1.2]. The effect of 
plant parameter variations may be incorporated into the 
unstructured uncertainty; however, the directional 
(phase) information associated with the structured 
uncertainty is lost and the result is an overly 
conservative design. Although these methodologies 

This research was supported by the NASA Ames and 
Langley Research Centers under grant NASA /NAG 2-297 , 
by the Office of Naval Research under ' contract 
GNR/N000M-82-K-0582 (NR 505-003), and by the 
National Science Foundation under grant 
NSF/ECS-82 10960. 


guarantee nominal closed-loop stability, 
stability-robustness with respect to the unstructured 
uncertainty, and nominal performance, the Issue of 
robust performance is ignored. 

Techniques which deal directly -with structured 
uncertainty have been developed. Horowitz [3-5] has 
proposed the so-called Horowitz templates which 
represent, at a particular frequency, the gain and phase 
changes associated with parameter variations as a 
region on a Nichols chart. A loop transfer function is 
derived from the templates which ensures closed-loop 
stability and a certain amount of performance over the 
possible range of parameters. This procedure is 
extremely tedious for more than a few frequency points 
and a great deal of judgment is required on the part' of 
the designer. 

Sideris and Safonov [6j approach the problem of 
structured uncertainty by examining a plant template in 
the complex plane. A series of transformations at each 
frequency is performed which maps the -irregularly 
shaped region in the complex plane onto the unit disk. 
The directional properties of the uncertainty are 
eliminated and the transformed problem is essentially a 
design with unstructured uncertainty. Nevanlirma-Pick 
interpolation is used to find the compensator which is 
then transformed back in such a way that it is a 
solution to the original problem. This approach may be 
promising, however it appears that a great deal of 
computation is required at each frequency point. 

The problem of robust stabilization in the presence 
of parameter uncertainty is addressed by Khargonekar 
and Tamenbaum [7], The authors do not deal with 
performance issues directly and simultaneous 
variations in the poles and zeros cannot be considered 
within the present framework. Doyle [8-10] has 
developed a new mathematical quantity y. the 
structured singular value, to handle structured 
uncertainty and the problems of robust stability and 
robust performance. While it may be a valuable 
analysis tool [1 1,12], a feedback synthesis methodology 
based on y is not yet available. 

It is the purpose of this paper to develop techniques 
for analyzing the stability and performance of a SISO 
feedback system which contains a plant model with 
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parameter uncertainty and unmode led dynamics. Based 
on the analysis, a robust control design method which 
uses concepts from LQG/LTR and H 00 theory is outlined. 
The compensator guarantees nominal closed-loop 
stability and stability-robustness with respect to 
structured and unstructured uncertainty. In addition, 
the H 00 norm of the worst-case sensitivity function 
over the operating band is minimized. Since this norm 
provides an upper bound on ail possible sensitivity 
functions, the closed-loop system satisfies the 
condition of robust performance. 

2. Problem Formulation 

Consider the feedback system in Figure 1 with plant 
model g(s) and compensator k(s). The physical plant to 
be controlled is single-input, single-output, linear, and 
time-invariant. The transfer function g t (s), which may 

be of infinite order, fully describes the true plant. The 
low-frequency behavior of g t (s) is captured in a 

parameterized n tri order transfer function g(s,£), where 
£ represents the parameter vector. The actual value of 
£ is £*. Unstructured uncertainty results from the 
fact that the n tf1 order transfer function g(s,£*) cannot 
completely describe the (possibly) infinite-dimensional 
plant g^(s). Using the multiplicative form of the 

modeling error as in [I], the frequency-domain 
description of the unstructured uncertainty is 
represented by the stable transfer function [1+A(s)l. 
The true plant can be described as follows. 

g t (s) = g(s,£*)[ 1 ♦ A(s) ] ( 1 ) 

It is assumed that a frequency-dependent magnitude 
bound on the unstructured uncertainty is available, but 
the phase of A(jw) is completely unknown. 

1 A(ju>)| 1 A 0 (u>) (2) 

In any practical situation £* is unknown; however, 
the true parameter values are bounded by the known set 
8 (i.e. £* e 8 ). For design purposes a nominal 

parameter vector £ is chosen from the set 8 . 
Structured uncertainty, which arises when £ and £ are 
not identical, can be described in the frequency-domain 
using the multiplicative model error representation as 
in Equation (I). 

g(s,£*) = g(s,£)[ I + 8 (s ,£*,£) ] (3) 

From (!) and (3), the true plant g*(s) can be represented 

in terms of the low-frequency model, the parametric 
uncertainty, and the unmcdeled dynamics. 


g t (s) = g(s,£)[ l + 8(s,£*,£) 11 1 + A(s) ] (4) 

In order to design the compensator k(s) and to 
analyze the feedback system containing the true plant, 
information about the nominal model g(s,£), the 
structured uncertainty S(s,£\ £). and the unmodeled 
dynamics A(s) is required. The nominal mode! and the 
magnitude bound on the unstructured uncertainty A 0 (o)) 

are known a priori. The exact value of 8(s, £*,£) is 
unknown, but the set containing all possible values may 
be constructed from 8 . Using Equation (3) and 
replacing £* with £. S(s,£,£) may be computed for a 
given £. 

S(s .e.e) = (g(3,£)/g(s,£)) - I for £, £ e 8 (5) 

The above equation defines the set of structured 
uncertainty U s . 

U s (s.£) = { 8(s,£,£) | £ £ 8 ) for £ f 8 (6) 

U s and the set of transfer functions A(s) satisfying 

Equation (2) lead to the definition of the set of possible 
plant transfer functions G(s) which contains the true 
plant g t (s). 

G(s) = { g(s) | S(s,££)eU s (s.fl) , | A(jo>) | < A 0 (u>) } (7) 
where g(s) = g(s,£)[ 1 + S(s.££)][ 1 ♦ A(s)]. 


The design method in section 4 will find a 
compensator k(s) such that the closed- loop system is 
stable for all g(s) e G(s). In addition to robust 
stabilization, the compensator must achieve closed-loop 
performance and performance-robustness. The 
sensitivity transfer function s(s) = [1+g(s)k(s)]~’ 
evaluated along the jco-axis governs the 
command-following and disturbance-rejection 
performance of the feedback system in Figure I. The 
performance objective of the design method is to 
minimize the H 00 norm of the worst-case sensitivity 
function over the frequency range where the commands 
and disturbances have energy (i.e. the operating band 
Q 0 ), subject to closed-loop stability. By definition, the 

magnitude of the worst-case sensitivity function is an 
upper bound on the magnitude of the sensitivity 
function for all g(s) e G(s) and for all o>. Thus the 
feedback system satisfies the condition of robust 
performance with respect to the worst-case 
sensitivity. 
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3. Analysis 

This section covers the stability and performance 
analysis of the feedback system in Figure l. The 
compensator k(s) has been designed based on a plant 
g(s) e G(s). The loop transfer function t(s)=g(s)k(s) and 
the sensitivity s(s) are the functions of interest. Let 
t(s,fi) = g(s,fi)k(s) be the nominal loop transfer function. 
Then, 

t(s) = t(s,e)[ 1 ♦ S(s, fi.fi) ][ l + A(s) 1 (8a) 

= t(s.£) + t(s,fi)8(s,fi.fi) + t(s.fi)[ i +S(s,ai)]A(s) (8b) 
3.1 Stability 

At the very least, the compensator k(s) ensures the 
stability of the nominal closed-loop system. From the 
Nyquist criterion the number of encirclements of the 
critical point in the Nyquist diagram of t(jwfi) is 
known. At this point it must be assumed that g(s,fi) 
has the same number of unstable poles for all fi e 0. 
Since the structured uncertainty does not change the 
number of unstable poles and A(s) is stable, the Nyquist 
criterion requires the number of encirclements of the 
critical point to remain unchanged for the loop transfer 
function t(jw). 

A stability-robustness condi tion may be derived by 
examining the loop transfer function t(jw) in the 
complex plane at a particular frequency w (Figure 2). 
From Equation (8b). there are two perturbation terms 
added to t( jw.fi) which may alter the stability of the 
nominal loop (i.e. change the number of encirclements 
of the critical point). The perturbation caused by the 
structured uncertainty is t(jw,fi)S(jw,fi,fi) and the 
second term. t(jw,fi)[l + S(jw,fi,fi)]A(jw), is a result of 
the combined effect of structured and unstructured 
uncertainty. This perturbation is represented as a 
circle in the complex plane because the phase of A(jw) 
is unknown. 

Before a stability-robustness test can be derived, it 
must be assumed that the structured uncertainty does 
not alter the stability of the nominal loop (i.e. the 
perturbation t(jw,fi)S( jw.fi.fi) does not change the 
number of critical point encirclements of the nominal 
loop t(jw,fi) ). Note that this is a very restrictive 
assumption; however, it is only used to derive a 
stabiiity-robustness condition for design purposes and 
it does not limit the classes of uncertainty which may 
be considered. Under this assumption the loop transfer 
function t(jw) will be closed-loop stable if the 
distance to the critical point ( i.e. j ln(jw)j ) is 
greater than zero over all frequencies. If this is the 
case it is impossible for the perturbations due to 


uncertainty to change the number of critical point 
encirclements. From Figure 2, the distance to the 
critical point d c of the loop t(jw) may be computed. 

The closed-loop system may be unstable if the critical 
point is located in the interior of the circular region in 
Figure 2 (i.e. d c is negative). 

d c (w.fi,fi) = | l + t(jw,fi)[l+S(jw.fi,fi)]|- 

| t(jw,fi)[l + 8(jw.fi.fi)]A(jw) j (9) 

where S(jw,fi,fi) e U s ( jw.fi). j A(jw) | < A 0 (w) 

To perform a worst-case stability analysis the 
minimum value of d c (w,fi J Q) must be found by searching 

over the set U s (jw.fi) and by replacing | A(jw)| by its 
bound A 0 (w). 

(^(w.fi) = min { | 1 + t(jw,Q)[l + 8( jw.fi.fi)! |- 

SeU s | t(jw,fi)[l + S(jw,fi.fi)]| A 0 (w) } (10) 

Theorem 1 : Sufficient Condition for Robust Stability 

i 

If 3j4w.fi) > 0 for all w, then the true closed-loop 
system is stable. 

Proof : From Figure 2 and the above assumption, 

the number of encirclements of the critical 
point cannot be altered by the uncertainty 
if the inequality in Theorem 1 is satisfied. 

Note the fundamental difference between structured 
and unstructured uncertainty. In the case of structured 
uncertainty, directional (phase) information is 
exploited. That is, only structured perturbations which 
dec r ease the distance to the critical point are of 
concern, and perturbations away from the critical point 
increase stability-robustness and can be ignored. No 
phase information is available for the unstructured 
uncertainty. Therefore it must be assumed that the 
unstructured perturbations are always in the direction 
of the critical point. 

3.2 Performance 

From the definition of the sensitivity function. 
s(s)=[l + t(s)]-', performance can be interpreted from the 
Nyquist diagram in terms of the distance to the critical 
point. That is, 

|s(jw)| = l/d c (w.fi.fi) * (n) 
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The worst-case sensitivity function magnitude is found 
as an immediate consequence of Equation (! !). 

|s(jo>)| = l/d c (a>i) (12) 

Define a scalar measure of performance £ as the H 00 
norm of the worst-case sensitivity over the operating 
band. 

£ = max |s(ju>)j (13) 

(OfQ 0 

The objective of the robust design method is to 
minimize £ (i.e. maximize performance), subject to the 
closed-loop stability condition in Theorem 1. Since 
|s(ju))| < | sCjoo) ] for all g(s) £ G(s) and for all <o, the 
requirement of robust performance is satisfied. 

4. Robust Design nethod 

A method is proposed to design a compensator which 
guarantees closed-loop stability and performance in the 
presence of uncertainty (i.e. robust stability and robust 
performance). The controller minimizes the H 00 norm of 
the worst-case sensitivity function over the operating 
band. 

The concept of Nyquist-shaoing is introduced as an 
integral part of the robust design process. While 
loop-shaping techniques are concerned with tailoring 
the magnitude of the loop transfer function, 
Nyquist-shaping requires magnitude and phase 
information to construct a target Nyquist diagram. The 
target satisfies all stability and performance 
conditions and represents the desired nominal loop 
transfer function. For the given plant, the robust 
design method produces a compensator which yields a 
nominal loop transfer function t(jw,fi) whose Nyquist 
diagram approaches the target to any required degree of 
accuracy. 

It is assumed that the low-frequency, nominal plant 
model g(s£), the parameter set © (and hence U s ), and 

the magnitude bound on the unstructured uncertainty 
A 0 (w) are known to the control system designer. Then, 

the first step in the robust design method consists of 
finding the target Nyquist diagram via the 
Nyquist-shaping procedure. The target curve cannot 
have an arbitrary shape in the complex plane. Stability 
and analytic function theory require the target Nyquist 
diagram to satisfy the following constraints. 

( 1 ) Nyquist stabili ty cr i ter i on 

(2) Bode's Integral Theorem [131 

(3) Theorem 1 


Subject to the above constraints, the target must 
minimize the H 00 norm of the worst-case sensitivity 
function over the operating band. The task of 
incorporating these constraints into a standard 
optimization problem is a subject of current research. 
While not all of the details of the Nyquist-shaping 
algorithm have been rigorously formalized, the basic 
idea is to start with a Nyquist diagram which 
corresponds to a loop transfer function which is known 
to be closed-loop stable with respect to the structured 
and unstructured uncertainty in the given plant model 
(i.e. satisfies Theorem I). This initial curve in the 
complex plane is continuously deformed in such a way 
that the above constraints are met and the H 00 norm of 
the worst-case sensitivity function over £ Q 0 is 

monotonicalty decreasing. The procedure terminates 
when acceptable performance has been achieved, or 
when it is no longer possible to improve worst-case 
performance without violating the stability-robustness 
condition in Theorem 1. 

Finding the initial curve for the Nyquist-shaping 
algorithm is relatively straight-forward. Since current 

methodologies can handle unstructured uncertainty, the 
idea is to transform the original problem (with 
structured and unstructured uncertainty) Into one with 
just unstructured uncertainty. This must be done in a 
conservative manner so that the resulting compensator 
guarantees a stable closed- loop system. 

A new multiplicative uncertainty may be defined by 
lumping the structured and unstructured uncertainty 
together. 

A'(s,ai) = s(s.ai) ♦ [ i *s(s.fliMs) ( 1 4) 

The plant model g(s) is represented as follows. 

g(s) = g(si)[ l + A‘(s,&£) ] (15) 

By ignoring the directional properties of the structured 
uncertainty a frequency-dependent magnitude bound on 
A'(s.fi,£) is computed at each frequency. 

A' 0 (o>,£) = max [|S(jo>.iL&)| + | 1 + S(jo>,&fD| A 0 (ou) } (16) 
6£U S 

| A'(jto,&S)| < A' 0 ((b,S) for £ 8, for all u> (17) 

With the nominal piant model g(s^) and the bound on 
the unstructured uncertainty A'o(u).S). LQG/LTR [14] or 
another suitable control design methodology is used to 
find a compensator k 0 (s) which guarantees a stable 
closed- loop system. The Nyquist diagram of the loop 
transfer function to(ju>)=g(jwi)k 0 (ju)) is the start. ng 
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point for the Nyquist-shaping procedure. Note that this 
is an overiy conservative control system design because 
the directional nature of the structured uncertainty has 
not been used. The Nyquist-shaping procedure results 
in a target Nyquist diagram which should remove this 
conservatism and improve performance. 

The target curve will typically correspond to a 
high-order system. A finite-dimensional, 

parameterized loop transfer function must be obtained 
from the target Nyquist diagram. This can be 
accomplished by a least-squares fit to the magnitude 
and phase data at specific frequency points. The 
stability of the finite-dimensional target should be 
checked by Theorem 1. If the robustness condition is 
not satisfied, a frequency-weighted least-squares 
procedure is used to improve the transfer function fit 
in the frequency range where Theorem l was violated. 
Alternatively, a higher order transfer function can be 
used to fit to the target Nyquist diagram. 

The finite-dimensional target loop is used in the 
Formal Loop Shaping LQG/LTR methodology to arrive at 
a compensator k(s). In [15], Stein and Athans outline 
the framework for using LQG/LTR to recover arbitrary 
(stable, minimum phase) target loop transfer functions. 
That is. the LQG/LTR methodology can be used to find a 
compensator k(s) such that the magnitude of the loop 
transfer function g(si)k(s) matches the magnitude of 
the target loop. Unfortunately, this application of 
LQG/LTR with Formal Loop Shaping requires the plant 
model g(s£) to be stable and minimum phase. Research 
is being conducted to remove this restriction. 

5. Implications For Adaptive Control 

A very active search for a robust adaptive control 
methodology is being conducted and a trend in the 
literature is developing. Many researchers now believe 
that a robust adaptive control system must consist of a 
robust system identification algorithm coupled with a 
robust control design method [16-20], This philosophy 
has been referred to as adaptive robust control, in 
contrast to robust adaptive control. Compensator 
redesign takes place infrequently compared to the 
system sample period and only when more accurate 
information about the system can be provided by the 
identification algorithm. While the robust design 
method presented here will be useful in its own right 
for fixed-parameter compensator design, the goal Is to 
develop an on-line algorithm as part of a practical (i.e. 
robust) adaptive control system. In addition, the design 
method will provide the initial guess for the adaptive 
compensator. 

Over the years a great deai of attention has been 
paid to the development of specific adaptive 
algorithms; however, very little consideration has been 


given to' an issue at the heart of the adaptive control 
problem: what are the performance benefits of adaptive 
control? In theory an adaptive control system provides 
better performance with respect to a fixed-parameter 
compensator because more information about the 
physical plant is incorporated into the design process 
(on-line). 

However, robust adaptive control systems rely upon 
some combination of external persistently exciting 
signals (to ensure good identification), slow sampling 
(to provide stability robustness with respect to 
unmodeled high-frequency dynamics. [21]). and extensive 
real-time computation (to provide safety nets which 
turn off the adaptation when it exhibits instability. 
[22]). These robustify ing measures degrade 
command-following and disturbance-rejection 
performance and tend to neutralize the anticipated 
benefits of an adaptive compensator. In light of these 
circumstances it is imperative that the decision to use 
adaptive control, for a real engineering application, is 
based upon a quantitative assessment of the costs and 
the benefits of an adaptive system. The robust design 
method proposed here produces the nonadaptive 
feedback system which minimizes the H 00 norm of the 
worst-case sensitivity function over the operating 
band. This system may serve as a performance 
benchmark to which an adaptive control system is 
compared. 


6. Conclusions 

A frequency-domain analysis of the stability and 
performance of a SISO feedback system with structured 
and unstructured uncertainty has been performed. The 
crucial analysis parameter is the distance to the 
critical point in the Nyquist diagram. Directional 
information (in the complex plane) associated with the 
structured uncertainty is exploited to reduce 
conservatism. A new method is outlined to design 
linear, time-invariant compensators for SISO plant 
models characterized by parameter uncertainty and 
unmodeled dynamics. The resulting feedback system 
minimizes the H 00 norm of the worst-case sensitivity 
function over the operating band and is guaranteed to be 
closed-loop stable. The concept of Nyquist-shaping was 
introduced as an integral part of the design process. 

It is conjectured that the robust design method can 
be used on-line in an adaptive control context. However 
before the decision to use adaptive control is made, the 
control designer must have a quantitative measure of 
the performance improvement over the ’best" 
nonadaptive system. The robust design method provides 
a benchmark for the performance comparison. 
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Figure 2. Representation of t(joo) 
in the complex plane. 
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